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The contributions to the parallel session on heavy flavour production 
are summarized. In more than 30 presentations new theoretical devel- 
opments and new experimental results from the Tevatron, HERA, LEP, 
the B factories and vN scattering were reported. 

1 Introduction 

Although traditionally centered around deep-inelastic scattering (DIS) and 
structure functions, this conference series has become one of the prominent 
places to discuss strong interaction physics in a wider sense. In the heavy 
flavour case in particular, the focus is on understanding in a common QCD 
framework the charm, beauty and top cross sections in various production en- 
vironments. 

The experimental results reported in this session were obtained in pp, pN, 
uN, ep, "fp, 77 and e^e" collisions. Only a personal and biased selection from 
them can be mentioned in this summary. Since the status of heavy flavour 
physics at HERA has been well covered in the introductory talk [1], the em- 
phasis is here on new results and on other machines. 

The theoretical results presented in this session capture various recent de- 
velopments in the area of inclusive charmonium, open-charm, and open-bottom 
production, as well as the heavy-flavour contributions to the structure functions 
[2, 3]. 



2 Charmonium (Theory) 

Since the discovery of the J/^ meson in 1974, heavy quarkonium has provided 
a useful laboratory for quantitative tests of QCD and, in particular, of the 
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interplay of pcrturbativc and nonpcrturbativc phenomena. The factorization 
formalism of nonrelativistic QCD (NRQCD) [4] provides a rigorous theoreti- 
cal framework for the description of heavy-quarkonium production and decay. 
This formalism implies a separation of short-distance coefficients, which can 
be calculated perturbatively as expansions in the strong-coupling constant a^, 
from long-distance matrix elements (MEs), which must be extracted from ex- 
periment. The relative importance of the latter can be estimated by means of 
velocity scaling rules; i.e., the MEs are predicted to scale with a definite power 
of the heavy-quark (Q) velocity v in the limit v <^ 1. In this way, the theo- 
retical predictions arc organized as double expansions in q,, and v. A crucial 
feature of this formalism is that it takes into account the complete structure of 
the QQ Fock space, which is spanned by the states n = with definite 

spin S, orbital angular momentum L, total angular momentum J, and color 
multiplicity a = 1,8. In particular, this formalism predicts the existence of 
color-octet (CO) processes in nature. This means that QQ pairs are produced 
at short distances in CO states and subsequently evolve into physical, color- 
singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the 
limit u — > 0, the traditional CS model (CSM) [5] is recovered. The greatest 
triumph of this formalism was that it was able to correctly describe [6] the 
cross section of pp ^ J/iIj + X at the Tevatron [7], which had turned out to be 
more than one order of magnitude in excess of the theoretical prediction based 
on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from 
severe conceptual problems indicating that it is incomplete. These include 
the presence of logarithmic infrared divergences in the 0{as) corrections to 
P-wave decays to light hadrons and in the relativistic corrections to S'-wave 
annihilation [8], and the lack of a general argument for its validity in higher 
orders of perturbation theory. In the case of "/p J/ip+X in photoproduction, 
where next-to-leading (NLO) corrections in the CSM happen to be well defined 
[9], the NLO prediction encompasses the rather precise HERA data [10] within 
the theoretical uncertainties. 

In order to convincingly establish the phenomenological significance of the 
CO processes, it is indispensable to identify them in other kinds of high-energy 
experiments as well, using the very MEs extracted from the Tevatron data [7]. 
Attempts to do so range from very successful to inconclusive for the time being. 
Successful examples include the comparison of theoretical predictions of 77 
J/ip + X [11] and pp J/ip + X [12] with data taken by DELPHI [13] at LEP2 
and by PHENIX [14] at RHIC, respectively. Comparison of HERA [15] data of 
ep — > J/ip -|- X in DIS with NRQCD predictions [16] are generally favourable 
with respect to both normalization and shape, except for the distribution in 
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inelasticity z, where a shape problem familiar from photoproduction [17] is 
recovered. As for the J/ijj polarization inpp ^ J/^ + X, measurements at the 
Tevatron [18] tend to undershoot the theoretical predictions for direct [19] and 
prompt [20] production in the bins of highest transverse momentum (pr)- 

While the /cT-factorization [21, 22] and hard-comover-scattering [23] ap- 
proaches manage to bring the CSM prediction of — > J/ip + X much closer 
to the Tevatron data [7], they do not cure the conceptual defects of the CSM; 
NRQCD is still needed for that. However, they reduce the relative importance 
of the CO contributions at leading order (LO). Since NLO corrections allow 
for an incoming parton to gain a kx kick by emitting another parton before 
entering the hard scattering, one expects A;T-factorization effects to become 
much less important at NLO. While the evaluation of NLO corrections is well 
understood in NRQCD, the extension of the fcr-factorization formalism beyond 
LO is an open problem. The list of open problems in connection with the kx- 
factorization formalism also includes the determination of the proper evolution 
equation (BFKL, CCFM, DGLAP, or mixtures thereof), the establishment of 
a rigorous operator definition of the unintegrated parton density functions, the 
restoration of gauge independence for off-shell external partons, the avoidance 
of ad-hoc rcgularization to keep propagators away from poles, which are some- 
times used in the literature [22] , and the recovery of the coUinear-parton- model 
results for asymptotically large values of pt in practical applications [22]. 

The color evaporation model [24] provides an operative procedure to esti- 
mate the production cross sections of heavy quarkonia on the basis of the ones 
of open heavy-quark pairs. It is intuitive and useful for qualitative studies, and 
it also significantly improves the description of the Tevatron data as compared 
to the CSM [25]. However, it does not account for the process-specific weights 
of the CS and CO contributions, but rather assumes a fixed ratio of 1 : 7. 

3 New Charmonium measurements 

The analysis of inelastic J/i/i production with HERA I data is essentially com- 
pleted. At this meeting, ZEUS presented preliminary cross section results 
for charmonium production in DIS [26] which confirm last year's HI findings, 
namely that NRQCD with sizeable CO contributions reproduces the data rea- 
sonably well, in particular at higher pT- However, the dependence on the 
inelasticity variable z remains at variance with the calculations also in the new 
ZEUS data. 

New experimental testing grounds may provide further clues to clarify the 
validity of the NRQCD factorization approach. In addition to the 77 collision 
results [13] (Sec. 2), data on charmonium production measured by HERA-B, 
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at the B factories in the e+e~ continuum, or at the Tevatron in the low pt 
regime are entering the stage. 

The HERA-B collaboration has published a first measurement of R^^ , the 
fraction of J/ip mesons produced via intermediate Xc states, which lies consid- 
erably below CSM expectations, but matches the NRQCD predictions well [27]. 
The total statistics of about 3000000 J/tp and 20000 Xc mesons produced on 
different nuclear targets allows them, for example, to measure the cross section 
dependence on isotope number A. In NRQCD this is not the same for J/tp 
and Xc, since CS and CO pre-resonant states contribute with different strength 
and interact differently with nuclear matter [28]. A measurement of the A 
dependence of the Xc cross section also forms an important ingredient for the 
interpretation of heavy ion collision data. 

A rather new testing field is charmonium production in the e+e" annihila- 
tion continuum below the T(4S') resonance [29]. The cross section was found 
to be more than an order of magnitude higher than expected, but the mo- 
mentum spectrum does not exhibit an enhancement at large momenta, which 
would signal strong CO contributions. According to the preliminary Belle re- 
sults, (82 ± 20)% of the J/tp mesons are accompanied by another pair of charm 
quarks. The J/ip recoil mass spectrum, measured now in a sample correspond- 
ing to 101 pb~^ and shown in Fig. 1, exhibits several clear cc resonances. From 
these data, rates for many exclusive (double) charmonium modes were be mea- 
sured, or earlier measurements confirmed with better precision; in addition all 
exclusive open charm ground states have been observed recoiling against J/ip 
mesons. 

New possibilities also open up in familiar environments. The Tevatron 
experiments are vigorously entering back on stage with new data from run II. 
By the time of the DIS conference, CDF, for example, wrote 150 pb~^ to tape 
and presented analyses making use of about 70 pb"-'^, to be compared to the run 
I data set corresponding to 110 pb~^. With their upgraded detector, trigger 
and data acquisition systems kinematic ranges and channels can be probed 
which have been inaccessible so far in highest energy pp collisions [30]. One 
example are lowered dimuon trigger thresholds which allow to trigger on J/^* 
mesons decaying at rest and thus to measure their production cross section 
over the entire momentum range down to px = 0. The existing measurements 
of run I start at > 5 GeV and cover less than 10% of the total cross 
section. This allows one to test the NRQCD description over a wider range; 
it provides increased sensitivity to those CO MEs (^5*0 and '^Pj) which also 
govern J/ ip production at HERA and should enable more stringent tests of their 
universality. However, a complete picture and a meaningful statement on the 
validity of the NRQCD factorization approach will only emerge once the NLO 
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Figure 1: Invariant mass of the system recoiling against J/^p mesons produced 
in e~^e^ annihilation below the Upsilon resonance, measured by Belle. 

corrections, which are shown to be important in the HERA photoproduction 
case, become available for hadroproduction as well. This requires a huge effort 
on the theoretical side, which should bear fruit in the not so far future. 

4 Open heavy flavours (Theory) 

The QCD-improved parton model implemented in the MS renormalization and 
factorization scheme and endowed with nonperturbative fragmentation func- 
tions (FFs), which proved itself so convincingly for light-meson inclusive pro- 
duction [31], also provides an ideal theoretical framework for a coherent global 
analysis of D- [32] and _B-meson data [33], provided that /i 3> toq, where 
H is the energy scale characteristic for the respective production process and 
Q — c,b. Then, at LO (NLO), the dominant logarithmic terms, of the form 
Q,n,ra+i jj^n j with 71= 1,2,..., are properly resummed to all orders by 

the DGLAP evolution, while power terms of the form (mg//Lt^)" are negligibly 
small and can be safely neglected. In this massless- quark scheme or zero-mass 
variable- flavour- number scheme (ZMVFNS), which is sometimes quite wrongly 
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referred to as NLL approxim,ation,^ the Q quark is treated as masslcss and ap- 
pears as an active parton in the incoming hadron or photon, having a nonper- 
turbative parton density function. The criterion /x » tuq is certainly satisfied 
for e~^e~ annihilation on the Z-boson resonance, and for the photo-, lepto-, and 
hadroproduction of D and B mesons with px ^ mg. Furthermore, the univer- 
sality of the FFs is guaranteed by the factorization theorem [34] , which entitles 
us to transfer information on how c and b quarks hadronize to D and B mesons, 
respectively, in a well-defined quantitative way from e~^e~ annihilation, where 
the measurements are usually most precise, to other kinds of experiments, such 
as photo-, lepto-, and hadroproduction. In Refs. [32, 33], the distributions in 
the scaled D- and S-meson energy x = 2E/,/s measured at LEPl were fitted 
at LO and NLO using, among others, the ansatz by Peterson et al. [35] for the 
c ^ D and 6 ^ B FF at the starting scale /io = 2mQ. In the D*± {B+/B°) 
case, the e parameter was found to be £c = 0.0851 and 0.116 [32] {sb = 0.0126 
and 0.0198 [33]) at LO and NLO, respectively. We emphasize that the value 
of £ carries no meaning by itself, but it depends on the underlying theory for 
the description of the fragmentation process, in particular, on the choice of the 
starting scale jUq, on whether the analysis is performed in LO or NLO, and on 
how the final-state collinear singularities are factorized in NLO. 

In the traditional massive- quark scheme or fixed-flavour-number scheme 
(FFNS), the Q quark is treated in the on-mass-shell renormalization scheme, 
as if it were a massive Icpton in triplicate, and it only appears in the final state, 
but not as an active parton inside the incoming hadron or photon. There are 
no collinear singularities associated with the outgoing Q-quark lines that need 
to be subtracted and absorbed into FFs. In fact, in this scheme, there is no 
conceptual necessity for FFs at all, but they are nevertheless introduced in an 
ad-hoc fashion in an attempt to match the D- and B-meson data. However, in 
the absence of a subtraction procedure, there is also no factorization theorem 
in operation to guarantee the universality of the FFs [34] . By the same token, 
such FFs are not subject to DGLAP evolution. Thus, there is no theoretical 
justification to expect, e.g., that a single value of the Peterson e parameter 
should be appropriate for different types of experiment or at different energy 
scales in the same type of experiment. In other words, the feasibility of global 
data analyses is questionable in this scheme. Moreover, this scheme breaks 
down for pT ^ mq because of would-be collinear singularities of the form 
as In (pr/mg), which are not resummed. However, it allows one to calculate a 
total cross section, which is infeasible in the ZMVFNS. 

^The nonlogarithmic corrections of relative order as are fully included, except for terms 
that axe suppressed by powers of I rriQ / fjp J . 
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The attempt to split the D- and S-mcson FFs into so-caUcd pcrturbative 
FFs (PFFs) and nonperturbative remainders is interesting in its own right. 
However, detailed analysis for £'*='=-meson FFs [36] revealed that such a pro- 
cedure leads to deficient results in practical applications. On the one hand, 
at NLO, the cross section da/dx of e"'"e~ annihilation becomes negative in the 
upper X range, at a; > 0.9, where the data is very precise, so that a low-quality 
fit is obtained imless this x range is excluded by hand [36, 37]. On the other 
hand, the LO and NLO predictions for other types of processes, such as pho- 
toproduction in ep scattering, significantly differ [36], which implies that the 
pcrturbative stability might be insufficient. 

The idea [33] of performing a coherent analysis of LEPl and Tevatron 
data of inclusive B-meson production was recently revived using an unconven- 
tional scheme named fixed-order next-to-leading-logarithm (FONLL) scheme, 
in which the traditional result in the FFNS and a suitably subtracted result 
in a ZMVFNS with PFFs are linearly combined [38]. However, some degree 
of arbitrariness is inherent to this procedure, as may be understood by notic- 
ing that the ZMVFNS term is weighted with an ad-hoc coefficient function 
of the form p^/ -|- 25m^) so as to effectuate its suppression in the low-pT 

range and that this term is evaluated at p'j. = -^Jp^ + ttiq while the FFNS 

term is evaluated at pT- Since the FONLL scheme interpolates between the 
FFNS and the ZMVFNS with PFFs, it also inherits some weaknesses of both 
schemes detailed above. In particular, the negativity of the NLO cross section 
of e+e^ D/B + X in the upper x range forces one to exclude the data 
points located there from the fit. In Ref. [38], this is achieved by resorting to 
what is called there the moments method, i.e., the large-a; region is manually 
faded out by selecting one particular low moment of the FF, namely the one 
corresponding to the average x value, thereby leaving the residual information 
encoded in the data unused. 

A rigorous theoretical framework that retains the full finite-mg effects while 
preserving the indispensible virtues of the factorization theorem, namely the 
universality and the DGLAP scaling violation of the FFs, is provided by the 
general-mass variable- flavour-number scheme (GMVFNS) [34, 39]. In a nut 
shell, this procedure consists in explicitly performing the toq limit of the 
FFNS result, comparing the outcome, term by term, with the ZMVFNS result 
in the MS scheme, and subtracting the difference terms from the FFNS result. 
Owing to the factorization theorem [34] , the hard-scattering cross sections thus 
obtained can then be convoluted with nonperturbative D- and S-meson FFs 
extracted from LEPl data using the MS scheme [32, 33]. This is consistent 
because the finite-mg terms omitted in Refs. [32, 33] are negligibly small, of 
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order m^/m^^. In fact, a rcanalysis of the LEPl data in the GMVFNS should 
yield FFs that agree with the FFs of Refs. [32, 33] within the errors. The 
GMVFNS was recently implemented for direct [40] and single-resolved [41] 77 
collisions. In the case of 77 —t D*^ + X at LEP2, the inclusion of finite- 
rric effects was found to reduce the cross section by approximately 20% (10%) 
at pt = ^TUc (3mc) [40], i.e., their magnitude is roughly m^/p^, as naively 
expected. By analogy, one expects the finite-mt, terms neglected in Rcf. [33] 
to have a moderate size, of order 20% (10%) at pr = 10 GeV (15 GeV). This 
is considerably smaller than the scale uncertainty and appears insignificant 
compared to the excess of the CDF data [42] over the traditional NLO analysis 
in the FFNS [43]. 

The GMVFNS scheme can be further improved by implementing the proper 
kinematic constraints on the threshold behaviour. In the case of the structure 
function , this was done in Ref. [2] by appropriately introducing the variable 
X = a; (1 + 4mg/Q^) (ACOT(x) prescription). This was shown to enforce the 
proper threshold behaviour, to simplify the implementation, and to yield robust 
predictions. 

5 Heavy quark heavy hadron trEinsition 

In a number of talks the current understanding of heavy quark fragmentation 

(or hadronization) was reviewed. This enters into every comparison of produc- 
tion cross sections with theory, to be discussed in the subsequent section. 

Above all, fragmentation is very well constrained by experiment. The most 
precise data are nowadays those on B meson production in e+e" annihilation 
at the Z resonance [44, 45] . The compilation of LEP and SLD measurements 
in Fig. 2 shows impressive agreement between the different analyses, which are 
based on exclusive {D* lepton), semileptonic or inclusive vertex methods to 
tag the b hadron and infer its momentum from the partially reconstructed final 
state. The figure displays the unfolded spectrum of the momentum fraction of 
the weakly decaying hadron, normalized to the beam energy. It reflects both 
the effects of the hard QCD radiation processes and of the soft transition from 
the quark to the bound hadronic state. The experiments have also studied 
how well the data can be described when different parameterizations are used 
to model the soft phase, and arrive at the same ranking [46], namely that the 
Lund and the Bowler shapes are preferred with respect to, e.g. the Peterson 
parameterization. Of course this conclusion is relevant only when a leading 
order (LO) parton shower approach as it is implemented in the JETSET Monte 
Carlo program, is used to model the perturbative phase. 

At HERA, the charm fragmentation function has been probed [47]. Since 
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Figure 2: Fractional momentum distribution of the weakly decaying b hadron, 
measured in e^e^ annihilation at the Z resonance. 

the initial state does not constrain the primordial quark momentum, the beam 
energy has been replaced by the jet energy to normalize the charm hadron's 
momentum fraction x. Hard gluon radiation thus affects the observable x 
spectrum less than in the e+e~ case. The measured distributions obtained by 
ZEUS, ARGUS and OPAL, which are of course not expected to agree in detail, 
exhibit nevertheless the same gross features, and fitting a parton shower Monte 
Carlo prediction to the ZEUS data yields a Peterson parameter consistent with 
e+e" findings. 

An enhancement at low x, attributed to gluon splitting, is observed only 
in the OPAL data, in the ARGUS and ZEUS case the energies are too low. 
Gluon splitting to cc and bb pairs has been explicitly constrained by studying 
special (3 and 4 jet) event topologies at LEP [45, 48]. For charm, a new 
and quite precise ALEPH measurement and a new world average of gcc = 
(3.01 ± 0.33)% has been presented at this conference [48]. This is somewhat 
above theoretical expectation, but consistent within errors. The situation is 
similar for the ten times smaller g^i, value [45]; altogether this shows that also 
higher order corrections are well under control. 

Another way to test the universality of the non-pcrturbative quark-hadron 
transition is to compare the fragmentation fractions [47], i.e. the relative rates 
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with which the different weakly decaying heavy hadron states are being pro- 
duced. ZEUS measured now, in addition to D* , and also the produc- 
tion of D"*" and Ac (to represent the charmed baryons) and have thus determined 
the complete set, so that the rates can now be normalized to the sum and no 
recurrence to MC for the total charm cross section has to be made anymore. Ta- 
ble 3 compares the ZEUS results with HI measurements, showing consistency 
of the HERA data among each other, and with e^e~ results, which confirms the 
universality of characteristics like the occurrence of strange mesons, baryons 
or vector mesons. These results further validate the practice to use the clean 
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Figure 3: Fragmentation fractions for charm hadrons. 

and precise measurements from e+e^ annihilation to model and interpret ep 
data. They give experimental support to a rigorous scheme to relate the QCD 
theory of quarks to observable single hadrons, which was sketched above and 
which provides a powerful means to probe parton dynamics. 



6 Open c, b and t production 

From an industrious community at HERA and other machines, cross section 
data on open heavy flavour production keep coming in and providing new 
challenges to the common QCD framework. 

The two-photon analyses at LEP2 are being finalized, bb production, after 
the L3 and OPAL observations, is now also seen by DELPHI [49] to be in excess 
of the massive NLO QCD prediction by more than a factor of 2. Direct and 
resolved photon processes contribute here with similar strength. The result 
appears to be in contrast to the two-photon data on charm production which 
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are reproduced by theory, within errors. An inspection of the pr spectrum 
(Fig. 4) reveals that data he rather at the upper side of the error range of the 
theoretical prediction, in particular for small transverse momenta. One should 




Figure 4: D* cross section in two photon interactions at LEP2, compared with 
fixed order (massive) QCD calculations. 

bear in mind that the bb data are dominated by final states with px rrib, 
populating a region where the factorization ansatz may be at the limit of its 
validity. 

Thanks to their upgraded trigger capabilities, CDF was able to present the 
first open charm cross section data from the Tevatron to this conference [30]. 
The results, presently extracted from a rather small amount of luminosity, are 
shown in Fig. 5 in comparison with resummed QCD calculations in the FONLL 
scheme.. The data tend to be above the central prediction, but come closer as 
Pt increases and thus follow the behaviour of the theoretical scale uncertainty. 

Several new results were reported from HERA. HI charm photoproduction 
cross sections have been updated using the full pre- upgrade luminosity [50]. 
They rely on events with a tagged final state positron and have therefore 
smaller statistics than the ZEUS measurements. They are found to be above 
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Figure 5: Charm meson production rate at the Tevatron, normalized to re- 
summed QCD predictions in the FONLL scheme. 



the massive fixed order QCD predictions, again particularly at low transverse 
momenta, see Fig. 6. The ZEUS data [51], also shown, have meanwhile a "ter- 
rible precision" , as it may be seen from the theoretical point of view. They are 
in principle sensitive to the charm density in the photon, in the massless pic- 
ture, however the calculations using different parton density parameterizations 
differ most at low px where the theoretical uncertainties are much larger than 
these differences. 

In DIS, a number of new results from their ongoing heavy quark production 
analyses were reported by CHORUS [52]. For example, their data on charm 
meson production in charged current neutrino nucleon scattering improve the 
experimental constraints on the strange content of the nucleon. - HI measured 
first differential jet cross sections [53], for dijet events produced in association 
with D* mesons. Comparing with Monte Carlo predictions using different 
schemes for the evolution of the parton cascade, agreement is observed within 
expectation. 

The power of jet measurements is impressively demonstrated by the ZEUS 
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Figure 6: Charm meson cross section, measured by HI and ZEUS at HERA, 
compared to massive (labeled "NLO") and massless (labeled "NLL") NLO 
QCD calculations. 

analysis of charm dijet angular distributions in photoproduction, which has 
recently been finalized [47]. The cross sections (Fig. 7) are measured as a 
function of cos9* = tanh(7yjet i — fljet2) which can be reconstructed from the 
two jet rapidities and corresponds to the polar angle of the outgoing parton 
in the partonic rest frame in a 2 — > 2 process. This is done separately for 
events with reconstructed x-y above or below 0.75, i.e. in regions where either 
direct or resolved photon contributions dominate. First, the dependence on the 
absolute value of cos 9* reflects the shape of the hard QCD matrix element; for 
low x-y the steep behaviour indicates the predominance of gluon exchange in 
resolved processes, as also observed in inclusive jet data. Associating the D* 
meson with one of the jet tags a final state charm quark. The data show that 
for low x-y the latter is found predominantly in the photon direction (negative 
cos 9*). This would not be expected in gluon gluon fusion (and is not seen 
in the direct, high dominated by photon gluon fusion), but matches 

well the predictions based on diagrams with a c quark originating from the 
hadronic photon structure in the initial state. Such processes are included in 
the fixed order massive QCD calculation, but only as HO correction to direct 
photoproduction. In the massless picture, this is a leading order process, and 
large NLO corrections are indeed found in the single inclusive case. One should 
also note the rather large uncertainty of the prediction in the massive scheme in 
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Figure 7: Angular distributions of charm tagged dijet events at HERA; dia- 
grams for charm excitation and boson gluon fusion. 

this region. This example shows how constraints in HERA kinematics and the 
quark hadron correlation in heavy meson production can be used to identify 
corners of phase space where improvements in the calculations are needed. 

b production at HERA is omitted in this write-up. Interesting discussions 
on the different cross section definitions of HI and ZEUS and on their respective 
shortcomings took place [54, 55, 56], but the overall picture changed with new 
HI photoproduction results appearing at the 2003 summer conferences [57], 
with better precision and presented in a way allowing direct comparison to 
ZEUS. They show consistency between old and new data, and between exper- 
iments in the observable range; altogether they tend to be above expectations 
based on NLO QCD in the massive scheme. 

In hadroproduction, HERA B has appeared on the scene with a b cross 
section measurement [58], which is in agreement with QCD predictions. 30 
times more statistics has been collected, so one awaits the measurement of 
differential cross sections. 

Tagged b jet cross sections are in principle less affected by theoretical uncer- 
tainties arising from coUinear radiation than B meson or b quark cross sections. 
The latest DO data from run 2 [59] are consistent with run 1 results and are 
compared with MC predictions in Fig. 8. 

The increased energy of the upgraded Tevatron implies an extra bonus of 
30% in the expected top production rate ~ and the possibility to test its energy 
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Figure 8: b tagged jet cross section at the Tevatron. 



dependence. The first top cross section [60] from run 2 was presented and found 
to be consistent with such expectations, within large errors still, but marking 
a first step towards QCD tests with differential top distributions. 



7 Conclusion 



As for heavy-quarkonium production, the NRQCD factorization approach [4] 
represents a rigorous theoretical framework that is renormalizable, devoid of 
infrared singularities, and predictive. It generalizes the CSM curing its defects. 
However, NRQCD predictions for charmonium suffer from sizeable normaliza- 
tion uncertainties at LO, and it is an important task for the future to evaluate 
NLO corrections, both in as and v"^. At LO, there is room for nonperturbative 
ingredients, such as kx factorization, hard-comover enhancement, and shape 
functions, which tend to reduce the relative importance of the CO processes 
and ameliorate the well-known shortages of NRQCD to describe the large-z 
behaviour of inelastic photoproduction at HERA and the charmonium polar- 
ization in hadroproduction at large px at the Tevatron. The necessity for 
such nonperturbative effects is expected to be reduced as NLO corrections are 
included. 

As for open heavy-flavour production, the notorious excess of measurements 
of inclusive D- and i3-meson production in various experiments over the respec- 
tive FFNS analyses implemented with phenomenological FFs not conceptually 
necessitated by factorization is expected to be reduced to an acceptable level 
by adopting the GMVFNS, in which the massive quark also appears in the 
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initial state and would-bc coUinear singularities are subtracted to match the 
MS calculation for pT ^ toq. Then, the factorization theorem is at work, 
guaranteeing the universality and DGLAP evolution of the (nonperturbative) 
FFs, and a global data analysis is feasible. 

From the experimental side, one notes that many new results continue to in- 
crease the variety of channels with which the QCD framework of heavy flavour 
production is confronted. The Tcvatron definitely presented itself as an upcom- 
ing charm, beauty and top factory. The gross picture for open heavy flavours 
is that theory undershoots, in particular at low pr, and that error bands in the 
predictions should be taken seriously, as they arc in general being exhausted 
once data appear. HERA analyses continue to explore new domains of precision 
and phase space; they hold the potential to provide further guidance to refine 
strategies in perturbativc QCD, in particular with the enhanced potential of 
the upgraded collider and detector coming to production mode this year [61]. 
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